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The  kinetic  effects  of  low  temperature  non-equilibrium  plasma  assisted  CH4  oxidation  on  the  extinction 
of  partially  premixed  methane  flames  was  studied  at  60  Torr  by  blending  2%  CH4  by  volume  into  the  oxi¬ 
dizer  stream  of  a  counterflow  system.  The  experiments  showed  that  non-equilibrium  plasma  can  dramat¬ 
ically  accelerate  the  CH4  oxidation  at  low  temperature.  The  rapid  CH4  oxidation  via  plasma  assisted 
combustion  resulted  in  fast  chemical  heat  release  and  extended  the  extinction  limits  significantly.  Fur¬ 
thermore,  experimental  results  showed  that  partial  fuel  mixing  in  the  oxidizer  stream  led  to  a  dramatic 
decrease  of  0  concentration  due  to  its  rapid  consumption  by  CH4  oxidation  at  low  temperature.  The  prod¬ 
ucts  of  plasma  assisted  CH4  oxidation  were  measured  using  the  Two-photon  Absorption  Laser-Induced 
Fluorescence  (TALIF)  method  (for  atomic  oxygen,  0),  Fourier  Transform  Infrared  (FTIR)  spectroscopy, 
and  Gas  Chromatography  (GC).  The  product  concentrations  were  used  to  validate  the  plasma  assisted 
combustion  kinetic  model.  The  comparisons  showed  the  kinetic  model  over-predicted  the  CO,  H20  and 
H2  concentrations  and  under-predicted  C02  concentration.  A  path  flux  analysis  showed  that  0  generated 
by  the  plasma  was  the  critical  species  for  extinction  enhancement.  In  addition,  the  results  showed  that  0 
was  produced  mainly  by  direct  electron  impact  dissociation  reactions  and  the  collisional  dissociation 
reactions  of  electronically  excited  molecules  with  02.  Moreover,  these  reactions  involving  electron 
impact  and  excited  species  collisional  dissociation  of  CH4  contributed  approximately  a  mole  fraction  of 
0.1  of  total  radical  production.  The  present  experiments  produced  quantitative  species  and  extinction 
data  of  low  temperature  plasma  assisted  combustion  to  constrain  the  uncertainties  in  plasma/flame 
kinetic  models. 
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1.  Introduction 

Recently,  plasma  assisted  combustion  has  drawn  attention  for 
its  potential  to  enhance  combustion  performance  in  gas  turbines, 
scramjets,  and  internal  combustion  engines.  Furthermore,  exten¬ 
sive  efforts  have  been  made  to  understand  the  kinetic  processes 
in  plasma-flame  interactions  [1-17].  However,  plasma  assisted 
combustion  involves  strong  coupling  effects  not  only  from  the  flow 
and  heat  interaction  but  also  from  the  simultaneous  production  of 
active  radicals,  excited  species,  ions/electrons,  and  other  interme¬ 
diate  species.  Consequently,  one  of  the  major  challenges  is  to 
decouple  the  complex  plasma-flame  interactions  so  that  the  ki¬ 
netic  contribution  of  individual  species  produced  by  plasma  can 
be  understood  at  the  elementary  reaction  level  [15,16].  As  such, 
it  is  of  great  importance  to  (i)  develop  well-defined  plasma- 
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combustion  systems  that  can  provide  quantitative,  fundamental 
experimental  data  to  understand  the  effects  of  plasma  generated 
species  and  (ii)  develop  and  validate  plasma-assisted  combustion 
kinetic  mechanisms. 

Many  experimental  studies  have  been  carried  out  to  understand 
the  role  of  plasma  generated  species  on  ignition,  flame  speed,  and 
flame  stabilization.  For  ignition  studies,  the  reduction  of  ignition 
delay  time  by  non-equilibrium  nanosecond  pulsed  discharges 
was  reported  in  Refs.  [11,18].  The  results  showed  that  the  ignition 
delay  times  could  be  reduced  by  about  an  order  of  magnitude  with 
the  action  of  plasma.  An  increase  of  atomic  oxygen  (0)  concentra¬ 
tion  was  observed  with  plasma  and  explained  as  a  major  factor  for 
the  decrease  of  ignition  delay  times.  Recently,  a  major  kinetic  path¬ 
way  of  O  production  with  plasma  has  been  revealed  through  the 
quantitative  measurements  of  O  and  excited  N2  concentrations 
[8,19],  which  has  advanced  the  current  understanding  of  plasma 
assisted  combustion  at  the  elementary  reaction  level. 

The  efforts  on  quantitative  understanding  of  the  enhancement 
mechanism  through  the  flame/plasma  interaction  has  been  made 
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by  several  research  groups  [12,15,16,20].  Ombrello  and  co-workers 
have  demonstrated  experimentally  how  to  isolate  singlet  oxygen 
(02(a1Ag))  and  ozone  (03)  from  other  plasma  related  species  and 
shown  that  both  (at  concentrations  of  several  thousand  ppm)  can 
enhance  flame  speeds  by  a  few  percent  [15,16].  The  experimental 
data  further  revealed  that  the  current  kinetic  model  for  C^^Ag) 
significantly  over-predicted  the  enhancement  for  hydrocarbon 
flames,  due  to  the  lack  of  understanding  of  collisional  quenching 
rates  of  C^^Ag)  with  hydrocarbon  species.  Enhancement  of  flame 
stability  by  a  non-equilibrium  plasma  discharge  or  an  electric  field 
in  partially  premixed  systems  also  have  been  observed  for  lifted 
flames  [13,14,21].  Additionally,  different  explanations  for  plasma 
assisted  combustion  have  been  given,  based  on  the  production  of 
H2  and  CO  [13],  the  formation  of  active  radicals  such  as  OH  and 
O  [14],  as  well  as  the  ionic  wind  [21  ].  Unfortunately,  in  spite  of  rig¬ 
orous  experimental  investigations,  the  kinetic  contribution  from 
the  plasma/flame  interaction  has  not  been  fully  characterized 
due  to  the  complexity  of  either  their  interaction  or  the  experiment 
itself. 

Recently,  a  nanosecond  pulsed  discharge  system  has  been  inte¬ 
grated  into  the  oxidizer  side  of  a  counterflow  flame  configuration, 
in  order  to  obtain  a  simplified  plasma/flame  interaction  [17].  The 
experimental  results  revealed  that  the  extinction  limit  of  diffusion 
flames  can  be  significantly  extended  by  applying  the  discharge  on 
the  oxidizer  side  to  produce  O.  However,  considering  a  practical 
application,  where  the  flow  condition  is  controlled  by  the  intensive 
turbulent  mixing,  partial/full  oxidation  of  the  premixture  with 
plasma  would  be  one  of  the  important  factors  that  dominate  the 
flame/plasma  interactions.  In  this  case,  it  is  important  to  investi¬ 
gate  how  the  chemical  reactivity  of  the  premixture  can  be  induced 
or  promoted  through  the  plasma/flame  interactions,  particularly  in 
the  low  temperature  regime  (<1000  K).  Moreover,  due  to  a  dearth 
of  quantitative  data,  the  current  plasma  combustion  kinetic  mech¬ 
anism  for  ignition  [8]  remains  untested  for  low  temperature  CH4 
oxidation.  Consequently,  it  is  essential  to  provide  quantitative 
measurements  of  species  not  only  to  validate  of  current  kinetic 
models  but  also  to  advance  the  fundamental  understanding  of  ki¬ 
netic  aspects  of  plasma/flame  interactions. 

Accordingly,  the  goal  of  the  present  study  was  to  investigate  the 
effect  of  plasma  assisted  CH4  oxidation  at  the  low  temperature  re¬ 
gime  on  the  extinction  limits  of  diffusion  flames  with  partially  pre¬ 
mixed  oxidizer  in  the  counterflow  configuration.  The  species 
concentrations  have  been  measured  with  TALIF  for  O  and  Fourier 
Transform  Infrared  (FTIR)  spectroscopy  and  Gas  Chromatography 
(GC)  for  stable  intermediate  species  in  order  to  characterize  the 
performance  of  nanosecond  pulsed  discharge  in  a  partially  pre¬ 
mixed  oxidizer  stream  and  to  validate  the  current  plasma-combus¬ 
tion  kinetic  model.  This  model  has  been  validated  further  against 
the  experimental  measurements  of  extinction  limits.  Path  flux 
analyses  have  also  been  conducted  to  identify  and  understand 
the  important  CH4  oxidation  pathways  in  plasma  assisted 
combustion. 


2.  Experimental  methods  and  kinetic  modeling 

2.1.  Partially  premixed  counterflow  flame  system  with  nano-second 
pulsed  discharge 

A  schematic  of  the  experimental  system  is  shown  in  Fig.  1.  The 
setup  consisted  of  a  pair  of  counterflow  burners  that  were  located 
in  a  low  pressure  chamber  with  a  volume  of  81  L.  The  inner  diam¬ 
eters  of  the  reactant  and  co-flow  nozzles  were  20  and  28  mm, 
respectively.  The  burner  separation  distance  of  the  opposed 
nozzles  were  maintained  to  be  equal  to  the  nozzle  diameter, 
20  mm,  to  avoid  the  perturbation  of  the  boundary  condition  from 


Ch/Ar/He/CFU 


Fig.  1.  Schematic  of  experimental  setup. 


thermal  expansion.  The  co-flowing  inert  (Ar)  curtain  provided  the 
isolation  of  the  flame  from  the  environment  and  burner  cooling 
as  well.  A  honeycomb  plate  was  placed  inside  the  reactant  nozzles 
to  ensure  that  the  flow  field  was  uniform.  The  counterflow  nozzles 
were  made  of  quartz,  and  two  parallel  bare  metal  electrodes  (a 
copper  alloy  that  is  corrosion  resistant  to  avoid  the  catalytic  effect 
of  copper)  were  located  inside  the  tubes  of  the  upper  nozzle.  The 
electrodes  were  15  mm  x  22  mm  x  1  mm  (thick)  and  were  sepa¬ 
rated  by  10  mm  (inset  in  Fig.  1).  The  electrode  edges  were 
smoothed  to  avoid  the  local  concentration  of  electric  field  and 
dielectric  breakdown.  The  front  ends  of  the  electrodes  were 
10  mm  away  from  the  exit  to  avoid  electrical  interaction  with 
the  thermocouples.  The  experimental  pressure  was  60Torr,  the 
diffusion  time  scale  (approximately  1  ms,  estimated  by  P/D,  where 
l  was  the  electrode  thickness,  and  D  the  diffusivity  (~10-3  m2/s)  of 
02  at  300  K  and  60  Torr)  was  much  shorter  than  the  flow  residence 
time  (approximately  10  ms).  Thus,  the  diffusion  was  sufficiently 
fast,  and  the  effect  of  the  electrodes  on  the  flow  field  inside  the 
tube  was  negligible.  To  ensure  that  the  stagnation  plane  was 
formed  in  the  middle  of  the  two  burner  exits,  the  momenta  of 
the  two  reactant  streams  were  matched.  During  the  experiments, 
He  and  Ar  were  used  as  dilution  gases  for  the  partially  premixed 
oxidizer  stream  (He  can  improve  the  uniformity  of  the  discharge 
[22],  but  it  was  difficult  to  maintain  the  flame  if  He  was  used  as 
the  only  dilution  gas  due  to  the  large  thermal  diffusivity  at  low 
pressure.  Thus  in  our  experiments,  a  He  and  Ar  mixture  was  used 
as  the  diluent).  The  composition  of  the  oxidizer  stream  was  fixed 
at  02/Ar/He/CH4  (0.26:0.32:0.4:0.02)  to  observe  the  effect  of  CH4 
doping  in  the  oxidizer.  The  fuel  stream  was  CH4  diluted  by  Ar  (fuel 
mole  fraction  varied  from  0.2  to  0.4).  A  nanosecond  pulse  generator 
(FID  FPG  30-50MC4)  was  used  to  generate  the  non-equilibrium 
plasma  to  activate  the  partially  premixed  oxidizer  stream.  The 
pulse  generator  is  capable  of  producing  repetitive  32  kV  pulses 
with  pulse  duration  of  12  ns,  full  width  at  half  maximum  (FWHM). 
The  frequency  of  the  pulse  generator  is  adjustable  from  1  to 
50  kHz.  A  typical  voltage  waveform  measured  by  a  LeCroy  high 
voltage  probe  (PPE20KV),  is  embedded  in  Fig.  1.  The  current 
through  the  electrodes  was  measured  with  a  Pearson  Coil  (Model 
6585).  The  pulse  energy  supplied  to  the  discharge  was  estimated 
from  the  time  integration  of  the  voltage  and  current  profiles  and 
found  to  be  approximately  1.27  mj/pulse.  The  characteristics  of 
voltage-current  were  independent  of  pulse  repetition  frequency. 
Therefore,  average  powers  were  calculated  to  be  5.1,  12.7,  25.4, 
38.1  and  50.8  W  at  the  respective  frequencies  of  4,  10,  20,  30  and 
40  kHz.  The  amplitude  of  the  voltage  was  fixed  at  5.3  kV  for  all 
tested  experimental  conditions,  as  shown  in  Fig.  1.  The  direct  im¬ 
age  of  the  nanosecond  repetitive  discharge  is  shown  in  Fig.  2a  at 
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Fig.  2.  Discharge  at  Ar/He/02/CH4  (0.32/0.4/0.26/0.02  by  volume),  d  =  20  mm,  /=  30  kHz,  P  =  60  Torr,  (a)  direct  image,  (b)  optical  emission  spectrum. 


a  pulse  repetition  frequency  of  30  kHz.  The  optical  emission 
spectrum  of  the  discharge  is  shown  in  Fig.  2b.  Figure  2b  shows 
strong  emissions  from  excited  O,  OH,  HCO,  CH,  He,  and  Ar,  which 
indicated  the  radical  formation  from  CH4  oxidation  during  the  dis¬ 
charge.  It  was  also  observed  that  the  strongest  emissions  were 
from  excited  Ar  and  O  (750  nm  and  777.4  nm,  respectively)  and 
the  intensity  of  all  emission  lines  increased  with  the  increase  of 
pulse  repetition  frequency.  Thus,  the  plasma  assisted  CH4 
oxidation  was  presumed  to  be  affected  by  the  pulse  repetition 
frequency,  and  the  excited  Ar  and  O  may  play  an  important  role 
in  the  CH4  oxidation.  The  intensity  of  excited  He  emission  was 
not  strong,  however;  this  is  because  the  excitation  energy  of  He 
is  relatively  high  (19.8  eV)  compared  to  that  of  Ar  (11.5  eV). 
Detailed  discussions  of  the  kinetic  process  will  be  given  in  the  next 
section. 

The  extinction  strain  rates  were  measured  by  fixing  the  fuel 
mole  fraction  and  increasing  the  flow  velocity  gradually  until 
extinction  occurred.  At  the  same  time,  temperatures  at  the  burner 
exits  were  measured  simultaneously  using  thermocouples.  One 
thermocouple  (380  pm  bead  diameter)  was  fixed  at  the  exit  of 
the  lower  fuel  nozzle,  and  another  movable  thermocouple  was 
used  to  monitor  the  temperature  at  the  exit  of  the  upper  nozzle. 
The  movable  thermocouple  was  coated  with  MgO  and  covered 
with  a  grounded  Nickel-Chrome  sheath  (with  1  mm  outer 
diameter)  to  remove  the  electromagnetic  interference  from  the 
discharge. 

The  major  products  (CH4,  CO,  C02,  CH20,  and  H20)  from  the 
plasma  assisted  CH4  oxidation  in  the  oxidizer  stream  were  mea¬ 
sured  by  using  an  FTIR  spectrometer  (Nicolet  Magna-IR  550).  A 
quartz  probe  with  an  inlet  diameter  1  mm  was  placed  axially  at 
the  oxidizer  side  nozzle  exit  and  was  attached  to  a  heated  line 
(to  avoid  H20  condensation)  of  the  FTIR  spectrometer.  The  pressure 
and  temperature  in  the  FTIR  system  were  held  constant  at  31  Torr 
and  393  K,  respectively.  For  H2  measurement,  a  micro  GC  system 
(INFICON  3000)  was  used.  The  calibrations  were  performed  by 
flowing  sample  gas  with  known  species  concentrations  through 
the  FTIR  or  GC  at  conditions  identical  to  those  of  the  experiments. 
The  relative  mole  fraction  uncertainties  of  sampling  measurements 
were  found  to  be  less  than  ±1%  for  CH4,  CO,  C02  and  ±5%  for  H20 
and  H2.  The  uncertainty  of  the  CH20  measurement  was  ±80  ppm. 
The  sampling  was  also  performed  by  using  a  quartz  probe  with 
6  mm  diameter.  The  results  were  identical  with  that  of  1  mm 
diameter  probe  which  indicated  the  uniformity  of  the  properties 
of  the  counterflow  system  in  the  center  region  and  the  validity  of 
one  dimensional  approximation. 


2.2.  Two-photon  Absorption  Laser  Induced  Fluorescence  (TALIF) 
system 


The  TALIF  method,  calibrated  with  Xenon,  was  used  to  measure 
the  absolute  O  concentration  at  the  oxidizer-side  nozzle  exit.  The 
TALIF  system  used  in  this  study  was  described  in  detail  in  our  pre¬ 
vious  study  [17],  and  thus  only  a  summary  of  this  method  is  pre¬ 
sented  here.  Ground  state  O  is  excited  by  absorbing  two  photons 
at  a  wavelength  of  225.7  nm.  The  transition  between  the  excited 
3p  3P  state  and  the  3s  3S  state  will  release  a  single  photon  at 
844.6  nm.  Xenon  can  be  excited  from  5p6  to  6p'  [3/2 ]2  with 
the  two  photons  at  224.31  nm;  de-excitation  to  6s'  [1/2]1  corre¬ 
sponds  to  fluorescence  at  834.91  nm.  The  number  density  of  O 
(N0)  was  calculated  using  the  following  equation  in  terms  of  the 
known  number  density  of  Xenon  (NXe)  [8], 


MXe) 

021  (O) 


f^(Xe)\ 

U(2)(0 )) 


1 

XW) 
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where  S0  and  SXe  are  the  observed  fluorescence  signals  for  O  and 
Xenon,  respectively,  a2i  =  the  fluorescence  quantum  yields 
(A2i  and  Q.  are  spontaneous  emission  and  quenching  rates,  respec¬ 
tively),  <j(2)  the  two  photon  absorption  cross  sections  of  Xe  and  O, 
F0(T)  the  atomic  oxygen  Boltzmann  factor  for  the  lower  level  of 
the  two  photon  absorption,  v,  the  photon  energies,  and  gND  the  neu¬ 
tral  density  filter  factor  (1.82  x  10“3  from  manufacturer).  These  val¬ 
ues  can  be  found  or  easily  calculated  based  on  the  data  in  Refs. 
[8,23]. 

An  Nd:YAG  laser  was  used  to  generate  532  nm  to  pump  a  tun¬ 
able  dye  laser  operating  at  ~573  nm.  This  573  nm  beam  was  fre¬ 
quency  doubled  and  mixed  with  the  1064  nm  beam  of  the 
Nd:YAG  laser  to  get  a  ~226  nm  beam  of  ~10pJ/pulse  required 
for  the  TALIF  diagnostics.  The  TALIF  signal  was  observed  using  an 
850  nm  bandpass  filter  of  40  nm  FWHM  and  a  Hamamatsu  photo¬ 
multiplier  (R636-10).  A  separate  vacuum  cell  was  used  to  obtain 
the  Xenon  atom  two-photon  fluorescence  spectrum.  Corrections 
for  the  transmission  of  the  vacuum  cell  and  bandpass  filter  were 
made  to  find  the  O  density  by  the  method  described  in  Refs. 
[8,23].  The  overall  relative  mole  fraction  uncertainty  of  the  TALIF 
measurement  was  approximately  ±40%. 


2.3.  Computational  approach  with  plasma-combustion  kinetic  model 

The  kinetic  model  used  in  this  work  was  an  integration  of  the 
air  plasma  model  [8]  together  with  USC  Mech  II  [24].  In  addition, 
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in  order  to  consider  species  dissociation  by  excited  Ar  and  He  and 
dissociation  of  CH4  by  electron  impact,  additional  elementary 
reactions  involving  the  dilution  gases  (He/Ar)  and  CH4  were  added 
(Table  1).  The  plasma  combustion  model  incorporates  a  set  of 
species  conservation  equations  for  number  densities  of  neutral, 
charged,  and  excited  species  produced  in  the  plasma,  0,  02,  03,  e, 
OJ,  0+,  0-,  02-,  02(a' Ag),  OzCb’X),  CX’D),  Ar(+),  Ar(3p5(2P3/2)4s), 
Ar(3p5(2P1/2)4s),  Ar(3p5(2P3/2)4p)  (Ar*  for  simplicity,  with  respec- 
tive  excitation  energies  of  11.5,  11.7,  and  12.9  eV),  He(+),  He(2s 
3S),  He(2s  ’S),  He(2p  3P),  He(2p  ’P),  He(3s  ’S),  and  He(4p  ’P)  (He* 
for  simplicity,  with  respective  excitation  energies  of  19.8,  20.6, 
21,  21.2,  22.9,  and  23.7  eV),  as  well  as  the  energy  conservation 


Table  1 

List  of  Ar/He/CH4  related  reactions  involved  in  Ar/He/02/CH4  mixture  discharge. 


Reaction 

Rate  (cm3  s  3) 

References 

e  +  CH4  =  e  +  CH3  +  H 

<7a 

[27] 

e  +  CH4  =  e  +  CH4(+) 

<7 

[28] 

e  +  CH4(+)  =  CH2  +  2H 

1.7  x  10“7  (300/T)0  5 

[28] 

e  +  CH4(+)  =  CH3  +  H 

1.7  x  10“7  (300/T)0  5 

[28] 

CH4(+)  +  02  =  CH4  +  02(+) 

5  x  10“10 

[28] 

Ar  +  e  =  Ar*  +  e 

o 

[28] 

Ar  +  e  =  Ar(+)  +  2e 

o 

[28] 

Ar*  +  02  =  Ar  +  20 

2  x  10-10 

[28] 

Ar(+)  +  02  =  Ar  +  02(+) 

1  x  10“10 

[28] 

Ar*  +  CH4  =  Ar  +  CH2  +  2H 

3.3  x  10“10 

[28] 

Ar*  +  CH4  =  Ar  +  CH  +  H2  +  H 

5.8  x  10“10 

[28] 

Ar*  +  CH4  =  Ar  +  CH3  +  H 

5.8  x  10-10 

[28] 

Ar*  +  CH4  =  Ar  +  CH2  +  H2 

5.8  x  lO”10 

[28] 

Ar(+)  +  CH4  =  Ar  +  CH3(+)  +  H 

6.5  x  10“10 

[28] 

Ar(+)  +  CH4  =  Ar  +  CH2(+)  +  H2 

1.4  x  10“10 

[28] 

He  +  e  =  He*  +  e 

<7 

[29] 

He  +  e  =  He(+)  +  2e 

<7 

[29] 

He*  +  02  =  He  +  02(+)  +  e 

1.5  x  10-11  T0-5 

[29] 

He*  +  O  =  He  +  0(+)  +  e 

1.5  x  io-11  T0-5 

[29] 

He(+)  +  02  =  He  +  0(+)  +  O 

0.6  x  10-11  T0’5 

[29] 

He(+)  +  03  =  He  +  02  +  0(+) 

0.6  x  io-11  T05 

[29] 

He(+)  +  02(f  Ag)  =  He  +  0(+)  +  O 

0.6  x  10-11  T05 

[29] 

He(+)  +  0(1D)  =  He  +  0(+) 

2.9  x  10-12T0'5 

[29] 

He  +  20  =  He*  +  02 

1  x  10“33 

[29] 

He  +  0(1D)  =  He  +  0 

1  x  10“10 

[30] 

He*  +  CH4  =  He  +  CH  +  H2  +  H 

5.6  x  10“13 

[31] 

a  The  rate  is  calculated  by  the  Boltzmann  solver  using  experimentally  measured 
cross-sections,  a. 


equation  for  prediction  of  the  temperature  of  the  mixture.  This 
set  of  equations  are  coupled  with  the  steady  state,  two-term 
expansion  Boltzmann  equation  for  the  electron  energy  distribution 
function  (EEDF)  of  the  plasma  electrons  using  experimentally 
measured  cross  sections  of  electron  impact  electronic  excitation, 
dissociation,  ionization,  and  dissociative  attachment  processes. 
The  Boltzmann  equation  calculates  the  rate  coefficients  of  the  elec¬ 
tron  impact  elementary  reactions  by  averaging  the  cross  sections 
over  the  EEDF.  The  kinetic  model  also  incorporates  chemical  reac¬ 
tions  of  excited  electronic  species,  electron-ion  recombination  and 
ion-ion  neutralization  processes,  ion-molecule  reactions,  and  elec¬ 
tron  attachment  and  detachment  processes.  Note  that  the  present 
model  does  not  solve  the  Poisson  equation  for  the  electric  field  and 
therefore  does  not  take  into  account  charge  separation  and  sheath 
formation  near  the  electrode.  So  in  the  present  work,  the  reduced 
electrical  field  ( E/N  ratio)  in  the  plasma  was  considered  to  be  an 
adjustable  parameter,  varied  until  the  CH4  concentration  (at  the 
nozzle  exit,  1  cm  away  from  the  end  of  the  discharge  region) 
matched  with  the  experimental  results  at  the  burner  exit.  CH4  con¬ 
centration  was  used  here  to  determine  the  E/N  ratio  because  CH4 
concentration  measurement  had  the  smallest  uncertainty  (1%  vs. 
40%  for  TALIF  of  O).  Therefore,  only  the  concentrations  of  CO, 
C02,  H20,  H2,  CH20,  and  O  from  the  simulations  were  compared 
with  the  experiments.  The  model  was  also  employed  to  calculate 
the  radical/atom  (OH  and  H)  concentrations  by  matching  the  O 
concentrations  from  the  simulations  with  the  corresponding 
experimental  values.  The  radical  concentrations  were  further  used 
as  the  boundary  conditions  for  extinction  limit  calculations  in  Sec¬ 
tion  3.2.  More  details  about  this  model  can  be  found  at  Refs. 
[25,26]. 


3.  Results  and  discussion 

3.1.  Effects  of  discharge  repetition  rate  on  species  concentrations  from 
plasma  assisted  CH4  oxidation 

FTIR  and  GC  instrumentation  were  used  to  quantify  the  major 
stable  products  of  the  low  temperature  plasma  assisted  CH4  oxida¬ 
tion  by  sampling  at  the  center  of  the  nozzle  exit  in  the  oxidizer 
stream  using  quartz  probe.  Typical  FTIR  absorption  spectra  are 
shown  in  Fig.  3a-f.  Figure  3a  shows  the  FTIR  absorption  spectra 


Wavenumber  (cm1)  Wavenumber  (cm1) 


Fig.  3.  FTIR  spectrum  with  different  pulse  repetition  frequencies  at  Ar/He/02/CH4  (0.32/0.4/0.26/0.02  by  volume),  P=  60  Torr;  (a)  no  plasma,  (b)/=  4  kHz,  (c )/=  10  kHz,  (d) 
/=  20  kHz,  (e)/=  30  kHz,  and  (f)/=  40  kHz. 
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without  the  discharge.  The  CH4  absorbance  spectrum  is  seen 
clearly  and  shows  no  effects  of  any  C02  and  H20  entrainment  from 
the  environment.  With  the  presence  of  the  plasma  discharge,  the 
absorbance  spectra  of  CO,  C02,  H20  and  CH20  appeared  due  to 
the  plasma  assisted  CH4  oxidation.  Moreover,  as  shown  in 
Fig.  3b-f,  with  the  increase  of  discharge  pulse  repetition  frequency 
(or  equivalently,  discharge  power),  the  concentrations  of  CO,  C02 
and  H20  increased  significantly,  and  at  the  same  time,  the  concen¬ 
tration  of  CH4  decreased.  The  dependence  of  the  measured  concen¬ 
trations  of  CO,  C02,  H2,  H20,  and  CH20  on  the  plasma  pulse 
repetition  rate  is  shown  in  Figs.  4-6.  In  order  to  validate  the  plasma 
flame  model,  the  predicted  species  concentrations  are  also  plotted 
in  Figs.  5  and  6.  In  all  the  FTIR  measurements,  it  was  confirmed  that 
the  carbon  closure  with  the  measured  species  was  95%,  demon¬ 
strating  that  the  measurements  captured  the  major  carbon  con¬ 
taining  species  in  the  plasma  assisted  CH4  oxidation.  It  was  also 
observed  that  carbon  black  deposited  on  the  electrodes  and  sam¬ 
pling  quartz  probe  during  the  experiments  and  that  the  minor  car¬ 
bon  loss  was  caused  by  the  formation  of  carbon  black  in  the  plasma 
discharge.  As  shown  in  Figs.  4  and  5,  with  the  increase  of  pulse  rep¬ 
etition  frequency,  the  CH4  concentration  decreased  significantly, 
and  nearly  a  mole  fraction  of  0.9  of  CH4  was  oxidized  at  pulse  rep¬ 
etition  frequency /=  40  kHz.  At  the  same  time,  the  concentrations 
of  CO,  C02  and  H20  increased  monotonically  with  /.  With  the  in¬ 
crease  of  CH4  oxidation  ratio,  the  oxidizer  temperature  at  the  noz¬ 
zle  exit  also  increased  significantly  (423  ±  3  K,  613  ±  5  K,  743  ±  5  K, 
843  ±  5  K  and  933  ±  8  K  at  /=  4  kHz,  10  kHz,  20  kHz,  30  kHz,  and 
40  kHz,  respectively).  Note  that  the  concentration  of  H2  after  low 


Pulse  repetition  frequency  (kHz) 

Fig.  4.  Dependence  of  CH4  concentrations  and  oxidization  ratios  on  pulse  repetition 
frequency. 


Pulse  repetition  frequency  (kHz) 

Fig.  5.  Dependence  of  species  concentrations  at  the  burner  exit  on  pulse  repetition 
frequency  (experiments:  solid  symbol  and  solid  line;  simulations:  open  symbol  and 
dashed  line). 


Pulse  repetition  frequency  (kHz) 

Fig.  6.  Dependence  of  species  concentrations  at  the  burner  exit  on  pulse  repetition 
frequency  (experiments:  solid  symbol  and  solid  line;  simulations:  open  symbol  and 
dashed  line). 


temperature  plasma  assisted  CH4  oxidation  remained  very  low 
(less  than  2000  ppm).  The  peak  concentration  of  H2  occurred  at  a 
pulse  repetition  frequency  /=  20  kHz.  The  low  value  of  peak  H2 
concentration  compared  to  H20  concentration  suggests  that  the 
plasma  assisted  combustion  was  not  dominated  by  fuel  reforming 
into  H2  and  CO,  as  will  be  discussed  later.  A  low  peak  concentration 
for  CH20  at/=  10  kHz  was  also  observed  during  the  experiments, 
as  shown  in  Fig.  6.  It  is  shown  in  both  Figs.  5  and  6  that  the  model 
captured  the  correct  trend  of  the  relationship  between  the  species 
concentrations  and  pulse  repetition  frequency,  but  the  kinetic 
model  over-predicted  CO  (up  to  44%),  CH20  (up  to  a  factor  of  6) 
and  H2  (up  to  a  factor  of  3),  and  under-predicted  C02  (up  to  89%) 
by  volume,  indicating  that  the  model  underestimated  the  overall 
reactivity  in  the  oxidizer  stream. 

In  order  to  identify  the  important  CH4  oxidation  pathways  and 
understand  the  kinetic  processes  in  the  plasma  assisted  CH4  oxida¬ 
tion  mechanism,  path  flux  analysis  was  performed  for  a  pulse  rep¬ 
etition  frequency  /=  40  kHz,  and  the  results  are  shown  in  Fig.  7. 
Methane  was  dissociated  to  CH3  mainly  by  H  abstraction  with 
OH,  O  and  H  produced  via  the  plasma  discharge.  Methane  also 
decomposed  to  CH3,  CH2  and  CH  via  collisions  with  an  electron, 
Ar(+)  and  Ar*  (reactions  (1  )-(6)).  The  formed  CH3,  CH2  and  CH  were 
oxidized  further  to  CH20,  HCO,  CO  and  finally  C02  as  shown  in 
Fig.  7a.  Among  those  reaction  paths,  a  total  of  a  mole  fraction  of 
0.12  of  CH4  was  decomposed  by  collisions  with  electrons,  Ar(+) 
and  Ar*.  Molecular  hydrogen  is  another  important  product  of  plas¬ 
ma  assisted  CH4  oxidation  because  it  can  affect  greatly  the  ignition 
and  extinction  processes  via  its  fast  oxidation  chemistry.  The  reac¬ 
tion  path  of  H2  is  shown  in  Fig.  7b.  The  dominant  formation  path  of 
H2  was  through  H  abstraction  reactions  with  CH4,  CH20,  H02,  HCO, 
CH2  and  CH3.  The  formation  of  H2  via  H  +  H02  =  H2  +  02  clearly 
indicates  the  nature  of  low  temperature  chemistry.  This  reaction 
plays  an  important  role  at  low  temperature  and  high  pressure 
when  H02  concentration  is  high  [32].  Another  formation  path  of 
H2  was  through  CH4  decomposition  by  Ar*,  but  only  1.2%  of  H2 
was  formed  through  this  reaction  path.  Moreover,  the  concentra¬ 
tion  of  H2  increased  at  low  pulse  repetition  frequency,  due  to  the 
increased  H  abstraction  reaction  by  H  radicals;  at  high  frequency, 
H2  concentration  then  decreased  via  increased  reactivity  at  higher 
temperature.  Therefore,  the  concentration  of  H  is  important  to 
determining  the  concentration  of  H2.  The  reaction  paths  of  H  and 
O  are  also  shown  in  Fig.  8a  and  b.  Only  11.7%  of  H  was  formed 
via  the  decomposition  of  CH4  by  electron  impact  and  excited 
species  collisions  through  reactions  (Rl),  (R2),  (R3),  (R4),  (R5). 
The  major  contribution  to  H  formation  was  from  reaction  CH3  +  O  = 
CH20  +  H.  Thus,  O  concentration  is  crucial  to  determining  H 
concentration. 
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Fig.  7.  Reaction  paths  flux  analysis  for  (a)  CH4  and  (b)  H2  at/=  40  kHz. 


Fig.  8.  Reaction  paths  flux  analysis  for  (a)  H  and  (b)  0  at/=  40  kHz. 


CH4  +  e  =  CH3  +  H  +  e 

(Rl) 

CH4  +  Ar(+)  =  CH3(+)  +  H  +  Ar 

(R2) 

CH4  +  Ar*  =  CH2  +  H  +  H  +  Ar 

(R3) 

CH4  +  Ar*  =  CH3  +  H  +  Ar 

(R4) 

CH4  +  Ar*  =  CH  +  H2  +  H  +  Ar 

(R5) 

CH4  +  Ar*  =  CH2  +  H2  +  Ar 

(R6) 

The  reaction  path  of  O  is  shown  in  Fig.  8b.  Approximately  54%  of 
O  was  formed  through  collisions  between  02  and  electrons/Ar*. 
Among  all  the  three  major  0  production  paths,  33.5%  was  from  di¬ 
rect  electron  impact  dissociation  of  02,  20.4%  from  02  collisions 
with  Ar*,  and  33.6%  from  the  H  +  02  =  OH  +  O  branching  reaction. 
Note  from  Fig.  8a  that  the  initial  H  radicals  for  the  branching  reac¬ 


tion  are  formed  from  the  reaction  of  O  with  CH3,  which  is  itself  pro¬ 
duced  by  the  direct  plasma  discharge  in  the  partially  premixed 
CH4/02  stream.  Therefore,  it  can  be  concluded  that  O  production 
by  plasma  is  the  major  source  to  initiate  the  low  temperature  plas¬ 
ma  assisted  CH4  oxidation.  Moreover,  the  partial  CH4/02  premixing 
increases  H  production  and  accelerates  the  chain  branching  pro¬ 
cess.  Once  O  is  generated  by  the  discharge,  it  is  consumed  rapidly 
by  CH4  and  its  intermediate  oxidized  and  dissociated  products 
(CH3,  CH2,  CH20,  HCO,  H02  and  H2).  It  is  shown  in  Fig.  8b  that 
the  major  consumption  path  of  O  is  to  generate  OH  and  H,  which 
are  important  to  the  CH4  oxidation  process.  OH  and  H  then  further 
react  with  fuel  and  its  fragments,  which  become  the  dominant 
reaction  paths  of  the  fuel  decomposition.  Only  a  small  amount 
(0.6%)  of  O  recombines  to  02  (Surface  recombination  reaction  of 
O  was  not  considered  in  this  study.  This  was  accounted  for  by  mea¬ 
suring  the  O  concentration  at  the  center  of  the  nozzle  and  compar¬ 
ing  with  the  simulations).  Note  that  the  reaction  path  from  O  to  03 
is  negligible  in  this  study  due  to  the  high  temperature  conditions. 
Furthermore,  the  concentrations  of  reactive  oxidant  ions  (Oj ,  0-) 
are  several  orders  of  magnitude  lower  than  O,  and  thus  their  effects 
are  also  negligible.  This  indicates  that  when  the  discharge  is  ap¬ 
plied  to  the  premixture  of  fuel  and  oxidizer,  O  can  be  consumed 
more  efficiently  by  the  fuel  and  its  intermediate  oxidation  and  dis¬ 
sociation  products  rather  than  through  recombination.  Based  on 
the  above  analysis,  it  is  clear  that  O  production  within  a  plasma 
discharge  in  a  partially  premixed  fuel  and  oxidizer  is  critical  for 
H2  and  H  formation  and  subsequent  CH4  oxidation.  Therefore,  in 
order  to  predict  other  species,  O  formation  through  the  collision- 
induced  02  dissociation  by  electrons  and  excited  species  needs  to 
be  accurately  predicted.  In  addition,  to  improve  the  CH4  oxidation 
efficiency  (the  consumption  efficiency  of  O),  the  plasma  discharge 
needs  to  be  implemented  in  a  partially  premixed  lean  fuel/oxidizer 
mixture. 

In  order  to  determine  the  concentration  of  O  and  compare  with 
the  kinetic  model  prediction,  the  Xenon  calibrated,  Two  Photon 
Absorption  Laser  Induced  Fluorescence  (TALIF)  method  was  used. 
The  relationship  between  O  concentration  and  pulse  repetition  fre¬ 
quency  (for  a  flow  velocity  of  200  cm/s)  with  and  without  fuel  pre¬ 
mixing  is  plotted  in  Fig.  6.  For  plasma  discharge  in  a  partial 
premixture  of  fuel  and  oxidizer,  it  is  shown  that  with  the  increase 
of  pulse  repetition  frequency,  the  O  concentration  increased,  but 
the  increase  was  very  low,  less  than  300  ppm  even  at  /=  40  kHz 
(note  that  the  TALIF  detection  limit  can  be  as  low  as  several  ppm 
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[8]).  The  low  rate  of  increase  of  O  production  versus /suggests  that 
the  cause  is  rapid  CH4  oxidation  at  high  plasma  repetition  fre¬ 
quency  by  the  induced  chain  reactions  (shown  in  Fig.  8b).  To  dem¬ 
onstrate  this  fact,  the  O  concentration  without  fuel  premixing  was 
also  measured  at  the  same  condition.  Figure  6  shows  that  the  con¬ 
centration  of  O  increased  sharply  from  530  ppm  to  2700  ppm 
when  the  pulse  repetition  frequency  increased  from  4  kHz  to 
40  kHz.  This  result  confirmed  that  the  increased  CH4  oxidation  at 
high  plasma  frequency  was  due  to  the  increased  production  of  O. 
The  simulation  results  of  O  and  CH20  concentrations  are  also 
shown  in  Fig.  6,  together  with  the  experimental  measurements  of 
CH20  concentration.  The  calculated  concentrations  of  O  and 
CH20  intersected  with  the  experimental  results,  but  the  deviations 
became  large  at/=  40  kHz. 

As  analyzed,  the  model  under-predicted  the  concentration  of 
C02,  and  over-predicted  CO  and  H2,  which  indicated  an  overall 
lower  reactivity  of  the  mixture.  For  O  and  CH20  prediction,  the  ki¬ 
netic  model  significantly  under-predicted  their  concentrations  at 
low  pulse  repetition  frequency  and  over-predicted  the  concentra¬ 
tions  at  high  pulse  repetition  frequency.  This  may  be  due  to  the 
missing  reaction  paths  involving  peroxide  species  at  low  tempera¬ 
ture  and  those  involving  interactions  of  excited  species  with 
hydrocarbon  species  and  its  products.  For  example,  the  H02  +  X 
(X  =  H,  OH)  reaction  rates  have  large  uncertainty  at  low  tempera¬ 
ture,  and  the  OH  +  O  +  M  =  H02  +  M  reaction,  which  is  not  included 
in  the  current  combustion  mechanism,  may  become  important  at 
low  temperature  [32].  In  addition,  a  fraction  of  O  atoms  was  ex¬ 
cited  to  0(1D),  which  is  much  more  reactive,  especially  at  low  tem¬ 
perature  conditions.  But,  unfortunately,  reactions  involving  0(1D) 
and  hydrocarbon  species  are  not  well  characterized  and  are  absent 
in  the  current  kinetic  model  (except  0(1D)  +  H2  =  H  +  OH  [33], 
which  is  107  faster  than  O  +  H2  =  H  +  OH  at  300  K).  The  formation 
of  CH20  was  primarily  from  the  reaction  CH3  +  O  =  CH20  +  H.  If 
the  concentration  of  O  was  over-predicted  (or  under-predicted), 
the  concentration  of  CH20  would  also  be  over-predicted  (or  un¬ 
der-predicted).  Likewise,  H20  can  react  with  0(1D)  via  the 
H20  +  0(1D)  =  20H  reaction  and  can  also  be  dissociated  by  elec¬ 
trons  and  Ar*  to  generate  OH  and  H.  The  generation  of  OH  can 
accelerate  the  oxidation  of  CO. 

The  kinetic  model  also  showed  the  formation  of  C2H6 
(~200  ppm),  but  the  C2H6  absorption  spectrum  was  not  observed 
in  the  FTIR  experiments.  Previous  research  in  a  lean  CH4/air  dis¬ 
charge  [34]  also  showed  no  C2H6  formation,  and  thus  the  reaction 
CH3  +  CH3  =  C2H6  was  over-estimated  or  the  consumption  path  of 
C2H6  was  underestimated  at  low  temperatures.  Furthermore,  it  is 
known  that  the  current  kinetic  mechanism  has  not  been  validated 
below  700  K.  More  quantitative  experimental  data  at  low  temper¬ 
ature  conditions  are  thus  required  to  improve  the  predictability  of 
the  kinetic  mechanism. 


In  this  experiment,  the  composition  of  the  partially  premixed 
oxidizer  stream  was  fixed  at  a  volume  ratio  of  0.02/0.26/0.32/0.4 
for  CH4/02/Ar/He.  The  fuel  stream  was  composed  of  CH4  and  Ar, 
and  the  CH4  mole  fractions  varied  from  0.2  to  0.4.  During  the  mea¬ 
surements,  the  fuel  mole  fraction  was  fixed,  and  the  flow  velocity 
of  both  sides  increased  gradually.  With  the  increase  of  flow  veloc¬ 
ity,  the  strain  rate  increased,  and  the  flame  had  less  residence  time 
to  complete  the  CH4  oxidation  and  thus  extinguished  when  the 
strain  rate  was  above  a  critical  value.  This  critical  strain  rate  at 
flame  extinction  was  recorded  as  the  global  extinction  strain  rate. 
Figure  9  shows  the  relationship  of  extinction  strain  rates  and  fuel 
mole  fractions  without  plasma  and  with  plasma  at  pulse  repetition 
frequencies  of  /=  4  and  10  kHz  (with  oxidizer  temperature 
T0  =  423  ±4K  and  613  ±  5  K,  the  oxidizer  side  flow  rates  were 
among  1.92-2.18  SLPM  and  2.22-2.41  SLPM  at  fuel  mole  fractions 
from  0.34  to  0.4,  respectively).  Extinction  strain  rates  without  CH4 
addition  (02/Ar/He  ratio  of  0.28/0.32/0.4)  were  also  presented  in 
Fig.  9  as  reference.  The  experiments  were  compared  further  with 
numerical  simulations.  The  simulations  were  computed  using  OPP- 
DIFF  of  the  CHEMKIN  package  [37]  with  a  modified  arc-length  con¬ 
tinuation  method  [38]  for  plug  flow.  All  simulations  were 
performed  by  setting  the  boundary  conditions  using  the  measured 
temperatures  and  species  concentrations  (using  FTIR,  GC  and  TALIF 
instruments).  The  concentrations  of  H  and  OH  were  estimated 
from  the  simulation  by  matching  the  computational  and  experi¬ 
mental  values  of  O  concentrations,  as  discussed  in  Section  3.1.  It 
was  noticed  that  the  extinction  strain  rate  was  sensitive  to  the 
CH4  mole  fraction  in  the  oxidizer  stream  because  this  determines 
the  heat  release  rate  in  the  reaction  zone.  Thus,  the  failure  to  up¬ 
date  the  species  concentrations  information  as  the  boundary  con¬ 
dition  may  cause  approximately  15%  larger  calculated  extinction 
strain  rates  in  the  simulations.  The  results  in  Fig.  9  show  that  with 
CH4  addition,  the  extinction  strain  rates  were  extended  signifi¬ 
cantly  and  the  oxidizer  temperature  after  the  plasma  activation 
was  higher  (613  I<  and  548  K  with  and  without  CH4  addition  at 
/=  10  kHz,  respectively).  It  is  also  shown  that  with  the  increase 
of  the  pulse  repetition  frequency ,/,  there  was  a  significant  increase 
of  extinction  limit  enhancement  and  the  measured  extinction 
strain  rates  agree  well  with  the  computed  values.  The  present  re¬ 
sults  indicate  that  the  plasma  generated  CH4  oxidation  plays  an 
important  role  for  the  extension  of  flame  extinction  limit.  This  is 
because  when  the  flame  approaches  extinction,  there  was  not  en¬ 
ough  residence  time  for  the  CH4  to  complete  the  oxidation  reac¬ 
tions  to  release  the  chemical  heat.  However,  plasma  can 
dramatically  accelerate  the  CH4  oxidation  to  release  the  chemical 
energy  at  low  temperature.  The  prompt  chemical  heat  release  ex¬ 
tended  the  extinction  limit. 

It  was  noticed  that  during  the  plasma  assisted  CH4  oxidation 
process,  part  of  the  CH4  was  reformed  into  H2  and  CO.  Molecular 


3.2.  Extinction  strain  rates  measurement  and  computation 


The  extinction  strain  rates  of  the  CH4/02  diffusion  flames  were 
measured  without  and  with  plasma  (at  /=  4,  10,  20,  30,  and 
40  kHz).  A  definition  of  global  strain  rate  [35] 
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was  adopted  and  compared  with  numerical  simulations.  Here  U  is 
the  speed  of  the  flow,  p  the  density  of  the  flow,  and  L  the  gap  dis¬ 
tance  between  the  two  nozzles.  The  subscripts  o  and  f  refer  to  the 
oxidizer  and  fuel  side,  respectively.  The  use  of  global  strain  rate 
was  found  to  be  appropriate  for  mechanism  validation  [36],  and 
the  global  strain  rates  were  used  in  this  study  for  both  experiments 
and  simulations. 
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Fig.  9.  Dependence  of  extinction  strain  rate  on  fuel  mole  fraction  at  the  fuel  side 
under  different  experimental  conditions. 
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hydrogen  can  extend  the  extinction  limit  through  its  fast  chemis¬ 
try,  and  previous  research  [13,20,39]  showed  that  plasma  induced 
fuel  reforming,  producing  a  H2/CO  mixture,  can  enhance  the  flame 
stability.  Therefore,  it  is  important  to  examine  how  the  fuel 
reforming  affects  the  extinction.  The  present  experiment  showed 
that  the  H2  formation  by  non-equilibrium  plasma  discharge  was 
low.  As  such,  the  plasma  reforming  effect  only  played  a  minor  role 
in  the  extension  of  the  extinction  limit.  Of  course,  the  H2  yield  de¬ 
pends  on  the  discharge  characteristics  and  fuel  concentration  in 
the  mixture.  It  is  necessary  to  understand  whether  fuel  reforming 
or  CH4  oxidation  is  more  effective  in  increasing  flame  stabilization. 

Because  plasma  can  cause  both  CH4  oxidation  and  reformation, 
simulations  of  the  individual  effects  of  plasma  assisted  CH4  oxida¬ 
tion  and  CH4  reforming  were  performed  and  compared  in  Fig.  10. 
Here,  part  of  the  fuel  (CH4)  in  the  oxidizer  stream  was  assumed 
to  be  ideally  oxidized  to  H20  and  C02  for  Case  1  and  reformed  to 
H2  and  CO  in  Case  2  under  a  constant  enthalpy  constraint.  It  is 
shown  that  extinction  was  enhanced  for  both  cases,  but  the 
enhancement  from  the  plasma  assisted  CH4  oxidation  (Case  1)  is 
much  larger  than  that  of  CH4  reforming.  This  is  because  when 
the  flame  was  approaching  extinction,  the  residence  time  was 
too  short  to  complete  the  reactions.  That  is,  CH4  and  its  reforma¬ 
tion  products  in  the  oxidizer  stream  cannot  react  completely  to  re¬ 
lease  the  chemical  enthalpy  in  the  reaction  zone  in  Case  2,  but  the 
chemical  enthalpy  was  already  released  in  Case  1  before  reaching 
the  reaction  zone.  Therefore,  plasma  assisted  CH4  oxidation  is  more 
effective  in  increasing  the  extinction  limit  than  CH4  reforming. 
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Oxidation  or  reforming  ratio 

Fig.  10.  Comparison  of  extinction  strain  rates  of  fuel  oxidation  and  fuel  reforming 
effect  under  a  constant  enthalpy  condition  for  the  oxidizer  stream  of  Ar/He/02/CH4 
(0.32/0.4/0.26/0.02  by  volume),  P  =  60  Torr. 


Fig.  11.  Dependence  of  extinction  strain  rate  on  pulse  repetition  frequency  at  a  fuel 
side  CH4  mol  fraction  of  Xf=  0.2,  an  oxidizer  side  volumetric  composition  of  Ar/He/ 
02/CH4  (0.32/0.4/0.26/0.02),  and  P  =  60  Torr. 


The  extinction  strain  rate  measurements  and  simulations  were 
also  conducted  at /=  20,  30  and  40  kHz  (with  oxidizer  temperature 
T0  =  743  ±  5  K,  843  ±  7  K  and  933  ±  8  K,  the  oxidizer  side  flow  rates 
were  1.66  SLPM,  1.89  SLPM  and  1.92  SLPM,  respectively,  and  with 
measured  and  calculated  boundary  conditions  as  described  above) 
by  fixing  the  fuel  mole  fraction  and  are  shown  in  Fig.  11.  With  the 
further  oxidation  of  CH4  (increased  f)  in  the  partially  premixed  oxi¬ 
dizer  stream,  the  extinction  strain  rates  increased  significantly. 
With  further  increase  of/,  the  experimental  and  simulation  results 
deviated,  which  may  indicate  that  there  are  additional  reaction 
paths  to  enhance  the  flame  extinction.  Nonetheless,  the  maximum 
deviation  was  small  at  approximately  10%.  Moreover,  the  peak 
concentration  of  H2  was  1751  ppm  at/=  20  kHz,  and  the  extinction 
limit  enhancement  from  H2  was  only  5.3%.  Thus,  the  dominant 
enhancement  mechanism  for  the  extinction  limit  is  plasma-in¬ 
duced  prompt  CH4  oxidation. 


4.  Conclusion 

The  present  experiments  employing  a  nanosecond  plasma  dis¬ 
charge  in  a  partially  premixed  counterflow  flame  showed  that  plas¬ 
ma  assisted  CH4  oxidation  at  low  temperature  (60  Torr)  can  extend 
the  extinction  limits  dramatically.  It  was  found  that  the  main  cause 
of  the  combustion  enhancement  was  due  to  a  kinetic-thermal  ef¬ 
fect  of  prompt  CH4  oxidation  and  heat  release  by  plasma  generated 
species.  Although  plasma  assisted  combustion  led  to  partial  CH4 
reforming  to  H2  and  CO,  numerical  simulations  showed  that  the  ef¬ 
fect  of  CH4  reforming  on  combustion  enhancement  was  mainly 
through  the  rapid  H2  chemistry  but  only  played  a  minor  role  com¬ 
pared  to  the  prompt  CH4  oxidation. 

The  results  also  showed  that  O  production  by  the  plasma  dis¬ 
charge  was  the  main  contributor  to  the  prompt  CH4  oxidation  at 
low  temperature  and  that  partial  premixing  of  fuel  in  the  oxidizer 
stream  led  to  a  dramatic  decrease  of  O  concentration  due  to  its  ra¬ 
pid  consumption  by  CH4  oxidation.  The  path  flux  analysis  showed 
that  O  was  generated  primarily  by  the  discharge  via  direct  electron 
impact  dissociation  of  02  and  02  collisions  with  electronically  ex¬ 
cited  species.  In  addition,  the  results  also  showed  that  the  electron 
impact  and  excited  species  collisional  dissociation  of  fuel  also  con¬ 
tributed  approximately  a  mole  fraction  of  0.1  of  the  radical 
production. 

The  extinction  strain  rates  from  experimental  measurements 
agreed  well  with  the  numerical  simulations  at  low  pressure.  Com¬ 
parison  of  measured  and  predicted  species  concentrations  in  plas¬ 
ma  assisted  CH4  oxidation  showed  that  the  low  temperature 
plasma  combustion  kinetic  model  reasonably  predicted  the  O 
and  CH20  formation,  but  the  model  over-predicted  the  concentra¬ 
tions  of  CO,  H20,  H2  and  under-predicted  the  concentration  of  C02. 
In  order  to  improve  the  predictability  of  the  plasma  flame  kinetic 
model,  the  elementary  reaction  rates  of  O  production  by  electron 
impact  and  excited  species  collision  need  to  be  revisited.  In  addi¬ 
tion,  more  experimental  validation  data  in  flow  reactors  of  CH4  oxi¬ 
dation  at  temperatures  between  300  K  and  700  K  with  radical 
addition  are  necessary. 
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